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We construct new solvable rational and trigonometric spin models with near-neigh- 
bors interactions by an extension of the Dunkl operator formalism. In the trigono- 
metric case we obtain a finite number of energy levels in the center of mass frame, 
while the rational models are shown to possess an equally spaced infinite algebraic 
spectrum. For the trigonometric and one of the rational models, the corresponding 



CN ■ eigenfunctions are explicitly computed. We also study the scalar reductions of the 

models, some of which had already appeared in the literature, and compute their 
algebraic eigenfunctions in closed form. In the rational cases, for which only partial 
results were available, we give concise expressions of the eigenfunctions in terms of 
generalized Laguerre and Jacobi polynomials. 
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Calogero-Sutherland (CS) models are one of the most extensively studied 
types of exactly solvable and integrable quantum Hamiltonians describing 
a system of iV particles in one dimension with long-range two-body inter- 
actions [1,2,3,4,5]. Apart from their intrinsic mathematical relevance [6,7,8], 
these models arise naturally in many different fields, such as Yang-Mills the- 
ories [9,10], quantum Hall liquids [11], random matrix theory [12,13,14,15], 
propagation of solitons [16], fractional statistics and anyons [17,18,19,20,21], 
etc. Several generalizations of CS models to particles with internal degrees of 
freedom ("spin") have been developed over the last decade using two main 
approaches, namely the supersymmetric formalism [22,23,24] and the Dunkl 
operator method [25,26,27,28,29,30]. Spin CS models are intimately connected 
with integrable spin chains with long-range position-dependent interactions, 
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like the celebrated Haldane-Shastry spin chain [31,32]. Indeed, when the cou- 
pling constant of a spin CS model tends to infinity the particles "freeze" in the 
classical equilibrium positions of the scalar part of the potential, thus giving 
rise to a spin chain of Haldane-Shastry (HS) type. This mechanism, usually 
called the "freezing trick" , was first used by Polychronakos [33] to construct 
the first integrals of the original HS spin chain by taking the large coupling con- 
stant limit of the corresponding integrals of the Sutherland spin model. These 
ideas have also been successfully applied to construct integrable and exactly 
solvable spin chains of HS type from other spin CS models [34,35,36,37]. 



In a recent paper [38], Jain and Khare proposed new solvable versions of the 
original Calogero and Sutherland scalar models featuring near-neighbors in- 
teractions. These new models are closely related to the so-called short-range 
Dyson models [39,40] in random matrix theory, in that the square of the 
ground-state wavefunction of the many-body system coincides with the joint 
probability distribution function for eigenvalues of the corresponding short- 
range Dyson model. Several generalizations of the Jain and Khare models ap- 
peared in a subsequent paper [41], including their extension to the BCn root 
system and to higher dimensions. The previous publications open a number 
of interesting questions, such as the existence of other solvable scalar models 
with near-neighbors interactions, or the construction of spin models with near- 
neighbors interactions and their corresponding spin chains of Haldane-Shastry 
type. These short-range spin chains with position-dependent interactions are 
of particular significance, since they occupy an intermediate position between 
the well-known Heisenberg chain (short-range, position-independent interac- 
tions) and the usual HS-type spin chains (long-range, position-dependent in- 
teractions). A first step in this direction is the recent work by Deguchi and 
Ghosh [42], in which several spin 1/2 models related to the scalar models of 
Jain and Khare were introduced and partially solved using the supersymmet- 
ric formalism. These authors also pointed out how to obtain the spin chains 
corresponding to these models by applying Polychronakos's freezing trick. 



In this letter we present three new families of solvable scalar and spin N- 
body models with near-neighbors interactions. By contrast with Ref . [42] , our 
approach is based on a modification of the Dunkl operator formalism and 
provides a wide range of totally explicit solutions for all values of the spin. 
The potentials of these spin models are given by 

Vi = 2a 2 cot (a* - Xi-i) cot(xi - x i+1 ) 
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+ 4a£-^^(a-S M+1 ), (lb) 

i \ X i X i+l) 

V 3 =u 2 r 2 + ^7 V + E? ~ — ^ (a - SWi) , ( lc ) 

^ (Xj — Xi-i){Xi — Xi + i) i [Xi — Xi + i) 

where r 2 = J2i% 2 and a,b > 1/2. Here and in what follows the sums run from 
1 to N, and we are identifying xjy + i with x 1 . The operators S^- permute the 
spin coordinates of the particles % and j. In other words, if \s±, . . . , s^) (with 
— M < Si < M, M being a half-integer) is an element of the basis of the spin 
space S, then 

Sij | ■ ■ ■ , Si, . . . , Sj, . . . } = | . . . , Sj, . . . , Si, . . . } . 

Note that can be expressed in terms of the fundamental SU(2M + 1) 
generators Sf, a = 1, . . . , 4M(M + 1), as = 1/(2M + 1) + £ a S^S?. 

The spin potentials (1) reduce to solvable scalar potentials by setting S^j+i 
to 1. In particular, the scalar reductions of the potentials (la) and (lc) are 
the models introduced by Jain and Khare. The potentials (la) and (lc) with 
M — 1/2 are similar to the spin 1/2 potentials introduced in Ref. [42], but 
differ from them by a spin- dependent term. The spin model (lb), as well as its 
scalar reduction, are both completely new. There is also a hyperbolic version 
of the trigonometric potential (la), obtained by replacing x k by ix k and V\ by 
-I 7 ,. 

The starting point in the solution of the models (1) is the introduction of the 
following second-order differential- difference operators, which play the same 
role as the quadratic combinations of Dunkl operators in the construction of 
CS models with spin [29,30]: 

Ti = E + 2a]T (z 2 di - z 2 +1 d i+ i) 

i i Zi — Z i+ i 

~ 2 «E r /^r , 2 (1 - (2a) 

T 2 = J2 Z i d i + 2a E ( Z i d i ~ z i+l9i+l) 

i i Z i~ Z i+l 

-°E ^ + z % (l-^i). ( 2b ) 
i {Zi — z i+1 ) 

T 3 = J2% + 2«E ft - 

i i z i~ Z i+1 

~ 2q E / \ v, (I-**,*!)- (2c) 

i \ Z i ~ Z i+l) 

Here di = d Zi , the operator Kij permutes the variables Zi and Zj, and we 
are again identifying z N+ i with z 1 . We shall also need in what follows the 
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first-order operators 



J- = 9 <> J = Y1 Z i d i ■ ( 3 ) 

i i 

If $ = Z)s/s( z )l s ), where z = (zi,...,zn) and s = (s 1? . . . , sn), is a state 
totally symmetric under permutations of both the spatial and spin coordinates 
of each particle, then = SV,$. Hence T e <3> = T e *$, e = 1,2,3, where 

T* = T £ \k ■ -ii->9 • The Hamiltonian H e = — J2idl + V £ of each of the 
models (1) can be obtained by applying a suitable gauge transformation and 
change of variables to a linear combination 

H e = cT* t +c^J_ + c J + E . (4) 

More precisely, we can write 

H e = fi-H e -pi' 1 , (5) 

Zi=C(Xi) 

where the constants c, c_, c , E , the gauge factor /i, and the function ^ are 
given in each case by 

1) c = 4, c_ = 0, c = 4(1 -2a), E = 2Na 2 , 

fi = Y[sm a ( Xl - x i+1 ), ((x) = e ±2ix , (6a) 

i 

2) c=-4, c_ = -2(26 + 1), c = 4c;, E = Nu{Aa + 26 + 1), 

y u = e-^ 2 n(^-x^ i rx J b , C(^) = x 2 , (6b) 

i 

3) c = -l, c_ = 0, c = 2cj, Eq = Nw(2a + 1), 

^ = e ~^ r " Il( x i ~ x i+iYi C(x) = x. (6c) 



The key idea in our approach to the solution of the models (1) is to find 
an increasing sequence of finite-dimensional linear spaces Tt® C 7i\ C 
invariant under the gauge Hamiltonian H t . Indeed, by Eq. (5) this implies that 
the corresponding Hamiltonian H e can be diagonalized in each of its invariant 
subspaces H™ = [iH^ _^ , n — 0, 1, The operators (2) and (3) preserve 

the space V n of polynomials in z of total degree at most n, for all non-negative 
integer values of n. In our recent work on spin CS models [29,30], the operators 
T cs analogous to T e also leave V n invariant, and in addition commute with the 
total symmetrizer under particle permutations A. This guarantees that the 
corresponding operators T* s preserve the space of completely symmetric spin 
functions A(V n <S> S). In the present case, however, the operators (2) do not 
commute with A, and hence it is not clear a priori whether the operators T* 
leave invariant any finite-dimensional space of spin functions. In fact, since 
J_ and J obviously commute with A, the results of Jain and Khare [38] for 
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the scalar case indicate that the operators Tj* and T 3 * possess at least "trivial" 
invariant spaces of the form Q<8>(A«S), where Q is a finite-dimensional subspace 
of the space of totally symmetric polynomials in z. As it turns out, each of 
the operators T* preserves a nontrivial subspace T e n C A(V n <g> <S) for all n, 
namely 

T\ = (f(r 1 ,T N -i,T N )A\s),g(T N - 1 ,T N )A(z 1 \s)),q(n,T N )A(zi ■ ■ -z N -i\s)) 

I fn = /jv-i,jv-i = gN-i,N-i = 911 — o) ) (7a) 
r 2™ = (/(Ti,7-2,r A r)A|s), 5 (r 1 ,r iV )A(z 1 |s)) | / 22 = = £atjv = o) , (7b) 
% n = (f(n, t 2 , T 3 )A\s),g{r u r 2 , r 3 )A(zi|s)), /i(n, r 2 )A(z?|s)), 

h{r 1 ,r 2 )A(z 1 z 2 \s')) | / 33 = ^33 = 0) . (7c) 

The proof of this statement, which is crucial for what follows, will appear in 
a separate publication. In Eq. (7), 

r k — ^ii - " " z i k 

ii<— <ifc 

is the /c-th elementary symmetric polynomial, /, g, h, and g are polynomials 
of total degree in z less than or equal to n, n — 1, n — 2, and n — iV + 1, 
respectively, and (for instance) = df/drk- The spin states \s) £ S are 
arbitrary, while \s') denotes a state such that the sum J2i \ s i,i+i) is totally 
symmetric, where |s^-) is defined by A^i^ls')) = Si<j ^i^'l 5 ^)- Note that 
the spaces (7) are also preserved by the operator J . Similarly, the operator 
J_ leaves the space invariant, while it preserves 7^ n and T£ provided that 
/jv-i = In = gN-i = gN = q = and = gN = 0, respectively. From 
Eqs. (4) jmd (6), it easily follows that the gauge Hamiltonians H e preserve the 
spaces given by 

n 1 — ii , n 2 — J-2 

By diagonalizing the Hamiltonian H t in its corresponding invariant spaces 
H™, n = 0, 1, ... , one can in principle construct an infinite sequence of exact 
eigenvalues and eigenfunctions, which shall be referred to as "algebraic" in 
what follows. We have found in this way all the algebraic eigenvalues of the 
spin models (1). We also present explicit expressions for the corresponding 
algebraic eigenfunctions, with the only exception of the spin eigenfunctions of 
the model (lc) not factorizing as /i/(r 1? r 2 , r 3 )A|s). In particular, we obtain all 
the algebraic eigenfunctions of the scalar reductions of the spin potentials (1), 
thus considerably extending the results of Refs. [38] and [41]. It should be 
noted, however, that the point spectrum could possibly include additional 
eigenvalues and eigenfunctions which are not algebraic. We shall now discuss 
in more detail each of the models (1). 



fn.T = On.T=0 ' 
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Case 1. The Hamiltonian Hi commutes with the total linear momentum P = 
— i J2k d Xk , so that it admits a basis of eigenfunctions with well-defined total 
momentum. In fact, since in this case t 1 n = exp ( ± 2HJ2k x k), multiplying 
an eigenfunction of Hi with energy E and total momentum p by t 1 n simply 
"boosts" its energy and total momentum. We shall take advantage of this 
fact to "normalize" the algebraic eigenfunctions to zero total momentum. It 
easily follows from Eq. (7a) that when this normalization is performed one 
obtains only a finite number of eigenfunctions of H 1 . In the scalar case, there 
are exactly four eigenfunctions with zero momentum, namely 



i IsinJ ^ 

Na 



^3 = H 



- + ^COS^ - Xj) 



2« + ^ i<3 



where X = J2i x % is the center of mass coordinate. Their respective energies 
are E (ground state), E ia = E + 4(2a -1 + 1/N), and E 3 = E + 8(2a + 1). 
These are essentially the solutions found in [43]. 

In the spin case, to each scalar eigenfunction (9) there correspond (^ M ^ N ^j 

factorized solutions of the form = ip n A\s), where \s) E S is arbitrary. 
There are three additional families of algebraic spin eigenfunctions with zero 
total momentum, given by 

=/^E{ COS ) (2( Xi - X))\ Si ) , 
% I sin ) 



where the spin states |sj) are defined by A^i|s)J = J2i z i\ s i), an d \s) £ S is 
any non-symmetric state. It should be noted that each of the spin states 
and has the same energy as the scalar eigenfunction tp n . 

Case 2. In this case the algebraic energies are given by E n = Eq + nco = 
Eq + Anuj, where the quantum number n — 0, 1, . . . is the degree in z of the 
corresponding eigenfunction of the gauged Hamiltonian H 2 . This follows easily 
from the fact that both T 2 * and J_ lower the degree, while J preserves it. As 
for Calogero's original model, the algebraic spectrum is equally spaced, but the 
spacing is twice the value suggested by the harmonic term. Unlike the usual CS 
models, the algebraic levels of this model have a well-defined thermodynamic 
limit, i.e., E n /N — > uj (4a +26+1) as — > 00. This property, which was already 
noted in Ref. [38] for the scalar reductions of the potentials (la) and (lc), is 
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in fact shared by all the models (1). In the scalar case, for each n > 2 there 
are two algebraic eigenfunctions with energy E n , namely 



n 



0,1, 



i 

+ n(n-l)aL«-\ur 2 ) 



n = 2,3, ... , 



where a = N(2a + 6+1/2) and is a generalized Laguerre polynomial of 
degree v. Note, in particular, that ip^ = fi has no nodes in the configuration 
space < Xi < ■ ■ ■ < xn, and is thus the ground state wavefunction. In 
the spin case, for each energy level E n we have first of all the factorized 
eigenfunctions of the form = ^>WA|s), k — 0, 1. In addition, for each 
n > 1 there is a family of genuine spin eigenfunctions 



#( 2 ) 



L«+_\{ur 2 )K(xl\s)) - ^-K-i("r 2 )A\s 



where \s) G S is non-symmetric. 



Case 3. This is probably the most interesting case, since the algebraic eigen- 
functions depend essentially on the three symmetric variables T\, r 2 , r 3 . As 
in the previous case, the algebraic energies are given by E n = Eq + hcq = 
E + 2nu, where n is again the degree in z of the corresponding eigenfunctions 
of H 3 . We shall begin, as usual, with the scalar case, for which we have been 
able to compute all the algebraic eigenfunctions in closed form. For each energy 
level E n with n > 3 there are two infinite families of algebraic eigenfunctions 
ipim, k = 0,1, with n = 21 + m. The first family is given by 



(o) 

Im 



HT™L- p (ur 2 )P^\t) , I, to = 0, 1, . . . 



J\ D (a,/J) 

m 
2 



(10) 



where t = 2Nr 2 /r 2 -l, a = N(a + 1/2) -3/2, (3 = (5{m) = l-m-N(a + l/2), 
p(a,P) j s a j acoD i polynomial of degree v, and [x] denotes the integer part of 
x. The second family reads 



(1) 

Im 



p, 



(q+3,/3) 



m — 3 i 
2 



WE 



+ 



2N 2 



¥>im(f) 



/, to — 3 = 0, 1, . . . , 
(11) 



where fi m is a polynomial of degree [to/2] given by 



to + 2a + 2 



m 



p^-2) (t) _ pi «+S^-l) (t) 



4a + 7 p (a+2,/J-l) 

±m -i 



m 



W + /frW. (12) 
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for even m, while for odd m we have 



+ m + 2 ' a t Pt + ^\t) - m + 2 " + 2 rfg »(t) . (13) 
m(m — 2) 2 m — 2 —r- 

Note that, as in the previous cases, ■^o°' ) = A* is the ground state. The above 
results show that for all N > 3 the scalar reduction of the Hamiltonian H 3 
possesses two different families of eigenfunctions of the form /iL^ /3 (wr 2 )p I/ (x), 
withp„ a homogeneous polynomial of degree v > 3, cf. Eqs. (10) and (11)— (13) . 
This was verified only up to v = 6 and N > v in Ref. [41]. 

In the spin case, the algebraic energies are of course the same as in the scalar 
case. As in the previous cases, each scalar algebraic eigenfunction gives rise 
to ( 2A ^" iV ) factorized spin eigenfunctions. The computation of the remaining 
algebraic spin eigenfunctions, which is considerably more involved than in 
Cases 1 and 2, is still in progress. Even without an explicit knowledge of the 
eigenfunctions, for each algebraic level E n one can easily compute the number 
of independent states of the form (7c) of degree n. If all the eigenfunctions of 
this model (and of its scalar reduction) were algebraic, the previous remark 
would imply that the degeneracy of the levels can be explicitly found. The 
method of Ref. [37] could then be used to compute the partition function of 
the associated spin chain 



where (£1, . . . , £tv) is an equilibrium of the scalar potential 



^ i ( x i ~ x i-l)( x i ~ x i+l) i ( x i ~ x i+l ) 



2 



This remark obviously applies to the rational model (lb) as well. Note, how- 
ever, that for this model one can also construct states of the associated spin 
chain by applying the freezing trick to the genuine spin eigenfunctions \l/^ 2 - ) 
presented above. 
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